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Abstract Semi-classical molecular dynamics simula- 
tions of small rare gas clusters in short laser pulses 
of 100 nm wavelength were performed. For comparison, 
the cluster response to 800 nm laser pulses was inves- 
tigated as well. The inner ionization dynamics of the 
multi-electron atoms inside the cluster was treated ex- 
plicitly. The simulation results underpin that at XUV 
wavelengths collisions play an important role in the en- 
ergy absorption. The generation of the surprisingly high 
charge states of Xe atoms inside clusters, as they were 
observed in the free-electron laser experiment at DESY, 
Hamburg, Germany [Wabnitz et ai, Nature 420, 482 
(2002)], is due to the reduced ionization potential of 
^ | atoms inside charged clusters, the ionization ignition 
1 mechanism, and collisions. 
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Clusters can absorb laser energy more efficiently than 
both gas targets and solids. This is because the clus- 
ter size is smaller than the skin depth so that all atoms 
experience the same laser field, despite the locally high 
particle density. Moreover, unlike solid targets there is 
no cold bulk that serves as a reservoir of cold electrons. 
Consequently, the laser-heated electrons cannot as easily 
escape from the cluster as they can from a solid surface. 
The electrons that absorb enough energy to leave the 
cluster as a whole contribute to the so-called "outer ion- 
ization." Other electrons are removed from their parent 
atom or ion ( "inner ionization" ) but remain trapped by 
the positive background of the cluster. These electrons 
and the ions form a nanoplasma that hydrodynamically 
expands or Coulomb-explodes. In experiments at wave- 
lengths > 248 nm highly energetic electrons |T], ions [21 
EE], photons 5 6,2], and neutrons originating from nu- 
clear fusion [H] were observed (see Refs. |5IT0*] for re- 
views). 

* E-mail: bauer@mbi-berlin.de 



Recently, Xe cluster experiments in the XUV regime 
have been performed at the DESY free-electron laser 
(FEL), Hamburg, Germany Surprisingly, for laser 
intensities of a few times 10 13 , 100 fs pulse duration, 
and 98 nm wavelength multiply charged Xe ions were 
observed in the laser-cluster experiments even for small 
clusters whereas isolated Xe atoms loose only one elec- 
tron under such conditions. 

The goal of this work is to shed some light on the 
mechanism underlying the generation of higher charge 
states inside clusters for laser parameters close to the 
FEL experiment. For that purpose molecular dynam- 
ics (MD) simulations have been performed. The inner- 
atomic dynamics is treated explicitly. The classical 
multi-electron atoms are rendered stable by introduc- 
ing a momentum-dependent potential. Mechanisms con- 
tributing to inner ionization such as field ionization, col- 
lisional ionization, and even non-sequential ionization 
and "shake-off" processes are all self-consistently incor- 
porated on a classical level. 

The paper is organized as follows. In Sec. |5] the nu- 
merical model is introduced. In Sec. OH the numerical re- 
sults are presented. Finally, we conclude in Sec. 0] 

Atomic units (a.u.) are used throughout unless noted 
otherwise. 



2 MD model of rare gas clusters 

The rare gas cluster consisting of N a atoms at the posi- 
tions Rj, 1 < i < N a , and Z "active" electrons per atom 
at the positions rj, 1 < j < ZN a in a laser field E(£) is 
modeled by the Hamiltonian 
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Apart from the usual terms describing the kinetic energy, 
the Coulomb interactions, and the interaction with the 
laser field in dipole approximation, there are the addi- 
tional potentials V aa and Vr accounting for the Lennard- 
Jones interaction between neutral atoms and the Heisen- 
berg uncertainty principle for the electrons inside the 
atoms, respectively. The purpose of introducing these 
potentials will become clear when we now describe how 
the initial, unperturbed cluster configuration is built up. 

(I) First, we seek an energetically optimal configura- 
tion for the neutral cluster without laser field. To that 
end we set Z — and 
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set to either unity or two in the chemical physics lit- 
erature. The potential minimum of V aa is at i? m i n = 
(2/£) 1//6 a and should be chosen close to the known near- 
est neighbor distance of the cluster under consideration. 
By varying the parameter D the overall strength of the 
Lennard- Jones potential V aa can be adjusted. 

Starting from a random atom distribution with the 
nearest neighbor distances > i? m ; n a local minimum in 
the energy landscape is obtained by propagating the 
atoms according the equations of motion 
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with a non- vanishing but small friction v and the nuclear 
mass M set to a small value so that relaxation occurs on 
an acceptable time scale. For small clusters (say, N a < 
10) it is quite likely that also the global minimum is found 
in this way as long as one starts with a reasonable guess 
for the cluster structure. Finding the global minimum 
for bigger clusters is far from trivial. Fortunately, the 
structures of Lennard- Jones clusters with N a < 150 are 
available in the literature |12II13| . 

(II) Secondly, we seek an electronic configuration 
for the single atom with Z active electrons and call 
this the "mother configuration." Since "classical atoms" 
with more than one electron are generally unstable, a 
momentum-dependent potential 
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is introduced [21 that enforces approximately the 
Heisenberg uncertainty relation when applied in the form 



V H (hj,Pij) = V(h 3 ,Pij,£,H,a H , Li ei ) 
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Here, f^ei — M/ (1 ~t~ -^0 ~ 1 is the reduced mass and Tij 
and pij are the absolute values of relative distance and 
momentum, respectively. 

The "hardness parameter" an governs how strictly 
the uncertainty relation rijpij > in is fulfilled. Big val- 
ues of an enforce it severely but also make the differen- 
tial equations of motion stiff (which is numerically un- 
favorable). We have chosen aj = 2. The parameter 
may be adjusted in such a way that the known essential 
ground state features of the system under study (e.g., 
the total energy or the ionization potentials) are prop- 
erly mimicked. 

Although we neglect the spin in the current study it 
is worth noticing that the Pauli principle could be mod- 
eled in a similar way so that an atomic shell structure is 
obtained (15116) . 

By solving the Hamilton equations of motion cor- 
responding to Q with a non-vanishing but small fric- 
tion v, the "mother configuration" of atoms or molecules 
with only a few electrons can be easily found. For "big- 
ger" atoms more advanced minimization routines have to 
be employed (151116) . Note that contrary to other classi- 
cal trajectory Monte Carlo (CTMC) methods where the 
classical atom is modeled by an ensemble of electrons 
moving on Kepler orbits, in our case the groundstate is 
stationary, that is r = p = but, owing to the Heisen- 
berg potential, p ^ 0. 

(Ill) Finally, the entire cluster is build up by tak- 
ing N a randomly rotated electronic mother configura- 
tions and attaching them to the naked ions. The ion 
positions are known from the first step. Since now all 
electrons "see" each other (as well as all the other ions) , 
another propagation with non- vanishing friction v is re- 
quired. During this step the electrons orient themselves 
in an energetically favorable way while the ions hardly 
move due to their huge mass M. Note that the cluster 
structure obtained in this way is not completely self- 
consistent because the ions are still sitting at the po- 
sitions determined in step (I) where just the Lennard- 
Jones potential was effective. In a fully self-consistent 
procedure the semi-classical cluster should be assembled 
without a Lennard- Jones potential at all. However, this 
is numerically much more demanding and may lead to 
an ion distribution which is physically less reasonable 
than the Lennard- Jones structure after step (I). 

We call the cluster configuration obtained after step 
(III) the "cluster mother configuration." In order to ob- 
tain meaningful results concerning the ionization dy- 
namics of the isolated atom or the cluster in a laser 
pulse, an ensemble of atoms and clusters had to be sim- 
ulated. The members of such an ensemble were con- 
structed by randomly rotating the cluster mother con- 
figuration. While this did not change the total energy 
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of the cluster or single atom, of course, the orientation 
with respect to the laser polarization axis changed, and, 
hence, so did the ionization dynamics. As expected, the 
ensemble-averaged entities converged more rapidly with 
respect to the time step for the integration of the equa- 
tions of motion (a fourth order Runge-Kutta scheme was 
used) as the individual particle trajectories did. Conver- 
gence of the latter may be accelerated by reducing an 
or introducing a "soft-core" smoothing parameter for 
the electron-electron and electron-ion interaction, that 
is, e.g., replacing \rj - | 1 by [(r,- - r z ) 2 + a^]" 1 / 2 
since all this loosens the requirement of small time steps 
during close encounters of the particles. 

3 Results 

We simulated small Xe clusters in intense laser fields. 
Each semi-classical Xe atom had Z — 3 active electrons 
that were supposed to mimic three of the six 5p electrons 
of Xe. Treating more than three electrons per atom is not 
necessary for the modest laser intensities we consider in 
this work. 

The other simulation parameters were chosen M = 
131 x 1836, D = 0.04, C = 1, #min = 8.0, a H = 2, 
£h = 1-75. The "mother configuration," i.e., a single Xe 
atom, consists of the three active electrons sitting on 
the corners of an equilateral triangle and the ion sitting 
in the center. The ion-electron distance was 1.6 for all 
three electrons. The total energy of the configuration 
was £i = —2.30 while the three ionization potentials 
were h = 0.32, I 2 = 0.80, and I 3 = 1.18. The first 
three ionization potentials for the real Xe atom are 0.45, 
0.77, and 1.18. The parameter = 1.75 was chosen 
to optimize I 2 and I3 since we are interested in multiple 
ionization and are particularly interested in the intensity 
regime where the removal of the first electron is given for 
granted for both the isolated atom and the atoms inside 
the cluster but higher charge states only occur in the 
cluster. Spin was neglected for if we also introduced a 
Pauli-blocking potential Vp of the form J5J into |Q we 
would have modeled the Li atom (with its closed s-shell 
plus a loosely bound valence electron) rather than three 
of the six 5p electrons of Xe. 

3.1 Isolated atom 

In Fig. Q] the charge state of the isolated Xe atom is 
shown vs the peak laser intensity of a T — 1764 (w 42 fs) 
laser pulse. The envelope of the electric field was linearly 
ramped up (and down) over T r = 331 (~ 8fs) while it 
was constant for the rest of the pulse. The two laser 
frequencies uji = 0.057 (corresponding to 800 nm) and 
u>h = 8a->; (corresponding to 100 nm) were used. 

As expected, for the same pulse duration and laser 
intensity the higher laser frequency ionizes the isolated 
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Fig. 1 Charge state of the isolated model Xe atom after 
a 42 fs laser pulse of different peak intensity. Low frequency 
case: u>i = 0.057, + symbols; high frequency case: uih = 8oj;, 
* symbols. Multiplication of the laser intensity by 3.5 x 10 16 
yields the laser intensity in the common units Wcm~ 2 . 

atom more efficiently than the lower one. For laser in- 
tensities I < 0.001 (« 3.5 x 10 13 Wcm~ 2 ) it is unlikely 
to observe charge states higher than the singly ionized 
Xe. This is in agreement with the results of Wabnitz et 
al. where at an intensity «2x 10 13 Wcm -2 for 12.7- 
eV photons (corresponding to 8.2 w;) and « 100 fs pulse 
durations also higher charge states than Xe + were found 
to be absent for isolated atoms. 



3.2 Small cluster 

In Fig. |21 the unperturbed Xe model cluster consist- 
ing of N a — 27 triply charged ions (black) plus the 
ZN a — 3 x 27 electrons (gray) is shown. The ion struc- 
ture is in agreement with Refs. [121113) . The electrons 
roughly maintain their triangular setup with respect to 
their parent ion. However, the triangles are in general 
not equilateral anymore due to the interaction with the 
neighboring atoms. 

The total energy of the unperturbed cluster amounts 
to £27 = —63.3 < N a £i = —62.1, i.e., the calculated 
cluster configuration is indeed energetically more favor- 
able than N a isolated atoms. The binding energy per 
atom is 0.04 which is a reasonable value for rare gas 
clusters. 

From Fig. ^ one can infer that for laser intensities 
/ less than 0.001 no charge states > 1 are expected in 
the case of isolated Xe atoms. Hence it is interesting to 
investigate whether higher charge states are observed for 
clusters. In Fig. [3] the results for the two laser frequen- 
cies 101 — 0.057, u>h — 8cJj and the two laser intensities 
J = 4x 10~ 4 (E = 0.02), I = 1.6 X 10~ 3 (E = 0.04) 
are collected for both the isolated atom and the cluster 
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Fig. 2 Unperturbed model Xe cluster configuration for 
N a = 27, Z — 3, and R m i„ = 8.0. Ions are drawn in black, 
electrons in gray. Three projections are also shown. 



consisting of N a = 27 atoms. The laser pulse duration 
and shape was the same as in Sec. 13.11 

The inner ionization of the cluster at the high fre- 
quency ujh (indicated by "inner (hf)" in the plots) is 
relatively high for both laser intensities. The inner ion- 
ization was calculated by counting the number of elec- 
trons that are more than 8 a.u. away from their initial 
position and therefore must have left their parent ion. If 
recombination plays an important role inner ionization 
is overestimated by this method. However, by looking 
at the particle trajectories we inferred that it is unlikely 
that an electron which left its parent ion is permanently 
trapped by another ion. 

Inner ionization is rather low for // and w\. This is 
because an important prerequisite is missing for these 
parameters, namely the efficient removal of the outer- 
most electron. It is clear that the very early inner ion- 
ization dynamics of clusters is the same as in the isolated 
atom case since there is not yet a strong influence of the 
other ions at that stage. In fact, the single atom ion- 
ization is small for Ii and oji (see Fig. dashed curve 
close to the bottom) so that the inner ionization of the 
cluster also remains quite modest. This is different for 
the higher laser intensity Ih (Fig. I^Jd) where the isolated 
atom looses the outermost electron. Consequently, in the 
case of a cluster a space charge builds up that yields fur- 
ther inner ionization ("ionization ignition," |17|). 

Outer ionization (calculated by counting the number 
of electrons which are farther away than r = 30) is al- 
ways smaller than the corresponding ionization of the 
isolated atom for the parameters in Fig. (see dotted 
curves). This means that a significant fraction of elec- 
trons that were removed from their parent ions do not 
make it to leave the cluster as a whole. This can be at- 
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Fig. 3 Number of removed electrons per ion for the two laser 
intensities (a) h = 4 X 10" 4 (E = 0.02) and (b) I h = 1.6 x 
10~ 3 (E = 0.04). The results for the isolated atom (dashed 
curves), the cluster outer ionization (dotted), and the cluster 
inner ionization (full) are plotted for the two frequencies ui = 
0.057 (If) and u h = (hf). 



tributed to the low quiver energies of the electrons inside 
the cluster so that no collective electron dynamics sets 
in that could help to absorb additional laser energy and 
to overcome the strong backholding field of the cluster 
ions ^Sl- Note that the ponderomotive energies for the 
parameters in Fig.[5]all lie between 4.8 x 10~ 4 and 0.12 
and are thus neither large compared to the photon en- 
ergy nor large compared to the binding energies. 

In Fig.0]the total energy (per atom) £(t) during the 
course of the laser pulse is plotted vs time. The laser 
and cluster parameters were the same as in Fig. |21 The 
net absorbed energy (per atom) \£(T) — £(0)| is for both 
frequencies and laser intensities greater for the cluster 
than for the isolated atom. As already mentioned in the 
discussion of the ionization degrees shown in Fig. |3 the 
absorbed energy is low for I\ and lo\ owing to the rather 
small probability for single ionization. For the higher in- 
tensity Ih instead, slightly more energy is absorbed at 
the low frequency u>i. This matches with the increased 
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Fig. 4 Total energy vs time for the cluster (full curves) and 
the isolated atom (broken curves). Laser and cluster param- 
eters were the same as in Fig. |3] 



outer ionization at low frequencies (cf. Fig. EJa) . More- 
over, the kinetic energy of the freed electrons depends 
on the laser frequency as well. 

In order to understand the absorption mechanism at 
short wavelengths it is worth to look at sample electron 
trajectories. In Fig. [S] one sees how electrons that at the 
end of the laser pulse are more than 30 a.u. away from 
the cluster center (and therefore contribute very likely 
to outer ionization) manage to escape from the cluster. 
As long as the distance from the initial position, r — r |, 
is less than ks 8 a.u., the electrons can be considered 
bound with respect to their parent ion. Zig-zag motion 
between distances from m 8 to ss 23 is due to collisions 
with other ions while electrons at distances > 23 left the 
entire cluster. 

For long wavelength and low intensity (a) one can see 
that the few electrons that made it to leave the cluster 
did so without collisions. This is because of the low prob- 
ability for the removal of even the outermost electrons: 
most of the atoms remain neutral so that it is unlikely 
that the few freed electrons encounter an ion. 




■■ 




















(O: 
h. to/ : 


20 


40 


6C 




Ir-rJ (a.u.) 



Ir-jjl (a.u.) 



Fig. 5 Sample trajectories of electrons which at the end of 
the laser pulse are more than 30 a.u. away from the cluster 
center. The absolute value of the distance to the electron's 
initial position, |r — ro|, is plotted for the four different cases 
(a) h,ui, (b) h,U)h, (c) Ih,wi, and (d) I h ,Uh- 



At higher frequency (b,d) it is seen from the erratic 
motion for 8 < |r — ro| < 23 that the electrons experi- 
ence several collisions with the cluster ions before they 
leave the cluster. Moreover, the electrons are temporar- 
ily trapped by other cluster ions. 

In plot (c), for twice the laser intensity and the low 
frequency, collisions are observed as well. This is because 
both the quiver velocity and the thermal velocity are 
small so that electron- ion collisions are likely (see, e.g., 
The outer electrons are removed efficiently at that 
laser intensity so that the ionization ignition mechanism 
starts (see the curve for inner ionization in Fig. Eb)- This 
indicates the eminent importance of the first freed elec- 
trons that trigger the entire cluster ionization process. 

From the trajectories we deduce that the efficient 
absorption of laser energy, as it is seen in Fig. is 
due to collisional heating or, in other words, inverse 
bremsstrahlung. Similar conclusions were drawn recently 
by other authors using different methods I2U1I2T1 . The 
increased inner ionization is due to (i) the reduced ion- 
ization potential because of the neighboring ions, (ii) the 
field generated inside the cluster by all other ions, and 
(iii) collisions. All three mechanisms are included self- 
consistently in our numerical treatment which, unfortu- 
nately, implies that it is difficult to separate the different 
contributions and to identify their relative importance. 

Although high charge states are generated inside the 
cluster during the course of the laser pulse it is not yet 
clear which charge states would arrive at the time-of- 
flight (TOF) detector in a real experiment. Upon the 
expansion of the cluster the electrons cool and some 
of them recombine, leading to a lower average charge 
state. Therefore we simulated the cluster disintegration 
as well. Since the cluster expansion happens on a pi- 
cosecond time scale we reduced the ion mass to M = 50 
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Fig. 6 Charge state distributions after the expansion of the 
cluster for uih and (a) (b) Ih- 



in order to make the simulation feasible [22] • While this 
mass reduction changes the time scale of the cluster ex- 
pansion it should not affect too strongly the hnal energy 
and charge state distribution as long as m/M <C 1. 

The result is shown in Fig. for the high frequency 
u>h and the two intensities Ii and Ih- In Fig. [Bp it is seen 
that despite the plasma cooling a significant amount of 
Xe 2+ ions survive the cluster expansion. In the case of 
the higher intensity Ih the doubly charged Xe ion is the 
most likely charge state and even triply ionized Xe atoms 
can be observed. 

Comparison of Fig. EJt with Fig. 1 in JT] (TOF spec- 
trum for N = 2-20) suggests reasonable agreement. 
However, a precise comparison with the experiment 1 1 1 j 
is not straightforward because of the not exactly known 
experimental laser intensity, the varying sensitivity of 
the TOF spectrometer with respect to different charge 
states, and focus effects. 



4 Conclusions 

Molecular dynamics was used to investigate small rare 
gas clusters in short XUV and 800 nm laser pulses. Inner 
ionization was treated explicitly. In order to render clas- 
sical, stable multi-electron atoms possible, a momentum- 
dependent potential that accounts for the Heisenberg un- 
certainty principle was introduced. Higher charge states 
in clusters than in the isolated atoms were observed, in 
accordance with experimental results. The mechanism 
underlying the efficient absorption of XUV laser energy 
was found to be inverse bremsstrahlung. 

The study of bigger clusters consisting of several hun- 
dreds or even thousands of atoms would be desirable. 
However, in order to keep the numerical simulations fea- 
sible one then has to give up the explicit treatment of 
the inncratomic dynamics. 
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